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The superconducting phase of hydrogen sulhde at Tc=200 K observed by Drozdov and collab¬ 
orators at pressures around 200 GPa is simple bee ImSm H3S, predicted beforehand by Duan et 
al, has experimental conhrmation. The various “extremes” that are involved - pressure, implying 
extreme reduction of volume, extremely high H phonon energy scale around 1500K, extremely high 
temperature for a superconductor - necessitates a close look at new issues raised by these char¬ 
acteristics in relation to high Tc. We use first principles methods to analyze the H3S electronic 
structure, particularly the van Hove singularities (vHs) and the effect of sulfur. Focusing on the two 
closely spaced vHs near the Fermi level that give rise to the impressively sharp peak in the density 
of states, the implications of strong coupling Migdal-Eliashberg theory are assessed. The electron 
spectral density smearing due to virtual phonon emission and absorption needs to be included ex¬ 
plicitly to obtain accurate theoretical predictions and current understanding. Means for increasing 
Tc in H3S-like materials are addressed. 


I. INTRODUCTION 

The recent discovery of superconducting hydrogen sul¬ 
fide under high pressure by Drozdov and collaborator 
and remarkably predicted a year earlier by Duan et al™ 
has reinvigorated the quest for room temperature su¬ 
perconductivity. The predicted structure has been con¬ 
firmed by x-ray diffraction studies W Shimizu that show 
that sulfur lies on a bee sublatticeP the protons cannot 
be seen in x-ray diffraction. The resistivity transitions 
were also confirmed by Shimizu. The experimental re¬ 
ports indicate critical temperatures up to Tc=203 K in 
the pressure range of 200 GPa, based on the resistiv¬ 
ity transition, the effect of magnetic field on Tc, on a H 
isotope shift of the right sign and roughly the expected 
magnitude,l2 and most recently the Meissner effect has 
been demonstrated.^^ 

In a success of predictive theory in this area, the mag¬ 
nitude of Tc in the 200 GPa pressure range was obtained 
from first principles calculation prior to experiment and 
confirmed by others,I^M^ so there can be little doubt that 
200 K superconductivity has been achieved in the struc¬ 
turally simple compound H 3 S, pictured in Fig. The 
finding that H vibrations provide the mechanism seems to 
confirm the suggestion of Ashcroft that dense hydrogen 
should superconduct at high temperatur^ however ev¬ 
idence is increasing that H-rich materials^ are substan¬ 
tially different and more promising than pure hydrogen 
until TPapressures can be reached. Early quantitative 
estimates^ of Tc for metallic H were in the 250K range; 
more recent valueJ^ in the range of several TPa lie in 
the 500-750K range. The phase diagram of this system 
is uncertain, however, due to the quantum nature of the 
proton. 

Although comprehensive calculations based on density 
functional theory (DFT) linear response formalism and 
Eliashberg theor}^^ have been reported and seem con¬ 
vincing, H 3 S turns out to be more intricate than the 
initial reports suggest. Using their self-consistent har¬ 



EIG. 1: Grystal structure of ImSm H 3 S. Nearest neighbor S-H 
bonds are shown. At this nearest neighbor level, the structure 
consists of two interleaved Re03 sublattices, displaced relative 
to one another by the body-centering vector (1,1,1)a/2. 


monic approximation, Errea et al. find substantial cor¬ 
rections due to anharmonicityPat 200 GPa, anharmonic- 
ity increases the characteristic frequency uiog by 3%, the 
electron-phonon interaction (EPI) strength A is decreased 
by 30% and the predicted value of Tc falls 22% from 
250 K to 194 K. Potentially important for further under¬ 
standing is their finding that anharmonicity shifts cou¬ 
pling strength to H-S bond stretch modes, from H-S bond 
bending (alternatively, H-H bond stretch) modes. 

Other basic questions have yet to be addressed. First, 
why are the electron-phonon matrix elements as large as 
they are? It is true that the main causative property be¬ 
hind the high Tc is the (understandably) high phonon fre¬ 
quencies that set the energy scale for Tc, but substantial 
electron-ion matrix elements are also required. Second, 
Flores-Livas et alt have investigated the energy depen¬ 
dence of the spectrum around the Fermi level, finding 


that it influences the theoretical predictions, which are 
overly optimistic when energy dependence is neglected. 
Both Akashi et al^ and Flores-Livas et al have solved 
the gap equations, providing theoretical predictions of 
the gap as well as Tc without using the Alien-Dynes equa¬ 
tion. This question posed by intricacies in the density of 
states (DOS) and the role of zero point vibrations has 
stimulated work by Bianconi and Jarlborg.E^ 

More fundamentally there is the question “why H 3 S? 
why sulfur?” Several H-rich materials have been stud¬ 
ied at high pressure (see references in Refs. HUE] and 
Bernstein et and although some are predicted to 

superconduct up to several tens of kelvins, H 3 S is a 
singular standout. Li et a/., for example, studied the 
H 2 S stoichiometry for stable compounds up to similar 
high pressures,^ finding a maximum Tc of “only” 80 K. 
There is little understanding so far of the microscopic 
cause of very high Tc, beyond the obvious expectation 
of higher phonon frequencies at high pressure; the origin 
of the large matrix elements remain obscure. Papacon- 
stantopoulos and collaborator^^ calculated the pressure 
dependence of matrix elements, finding increasing H scat¬ 
tering with increasing pressure. More basically one can 
ask, is there something special about sulfur, and the un¬ 
derlying electronic structure, that provides the platform 
for such high Tc? 

It is the last of these questions we address initially in 
this paper. An obvious feature for study is the strikingly 
sharp peak in the density of states N{E) due to two van 
Hove singularities (vHs) separated by 300 meV very near 
the Fermi level Ep. There is a large literature on the 
connection between peaks in N{E) and high Tc in the 
A15 class of materials^ a nd late r in the high temperature 
superconducting cnpratesj ^^ l ^^ * but their importance for 
H 3 S is unclear, van Hove singularities near the Fermi 
level can enhance N{Ef) and thus the EPC strength A 
due to increased number of available states to participate, 
but there are additional questions to address. 

The paper is organized as follows. Methods are de¬ 
scribed in a brief Sec. H. In Sec. HI the general electronic 
structure and the charge density near Ep are presented 
and discussed. Based on a Wannier function represen¬ 
tation of the bands, a minimal tight binding model is 
presented, with the intention of identifying the impor¬ 
tant features of the bonding and especially the deviation 
from free electron like density of states over much of the 
valence band. The two van-Hove singularities are identi¬ 
fied, quantified, and analyzed, and the relation between 
them is identified. In Section IV we address the peak 
in N(E) in the light of strong electron-phonon coupling, 
high frequencies, and thermal smearing. Sec. V presents 
scenarios for further increase in coupling strength A, and 
raising of Tc toward room temperature, in this and sim¬ 
ilar systems. A short Summary is provided in Sec. VI. 


II. METHODS 


Density functional calculations have been carried out 
using both the linearized augmented plane wave (LAPW) 
method based WIEN2k cod^^ and the linear combina¬ 
tion of atomic orbitals based EPLO code.l^ The PBE 
implement at ioiP^ of the generalized gradient approxima¬ 
tion (GGA) is used as the exchange correlation func¬ 
tional. The crystal structure of H 3 S is ImSm with a 
lattice constant of 5.6 a.u. corresponding to a pressure 
of 210 GPa. 




FIG. 2: Full density of states for WIEN2k (top panel) and 
FPLO (bottom panel) calculations. The bandwidth is ex¬ 
tremely large due to the broadening caused by short inter¬ 
atomic distances. The region of the van Hove singularities 
differ in detail. N{Ef) from FPLO is 0.44 states/eV, while 
from WIEN2k it is 0.63 states/eV. 


To study the vHs points near the Fermi level, a very 
fine k-mesh containing 8094 (8112) points for WIEN2k 
(FPLO) in the irreducible Brillouin zone is used. Sphere 
radii R for H and S are 0.97 uq and 1.81 uq respectively, 
with basis set cutoff determined by Rn^max = 6 . The 
results we discuss are insensitive to these choices. The 
tight-binding parameters we present were obtained as a 
two-center Slater-Koster simplification of a more exten¬ 
sive representation in terms of symmetry-adapted Wan¬ 
nier functions as implemented in FPLO. 
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III. ELECTRONIC STRUCTURE AND 
BONDING 


B. Minimal tight-binding model 


A. Electronic structure at 210 GPa 

The calculated lattice constaniPat 210 GPa, a=5.6 a.u. 
corresponding to a volume 58% of the zero pressure vol¬ 
ume in the same structure, is used in all calculations. The 
DOS N{E) on a broad scale from the two all-electron, 
full potential codes are presented in Fig. 2 . Because fine 
structure is of interest here, we compare a few results 
from two all-electron, full potential codes. The occupied 
bandwidth is 26 eV. Over the lowest 20 eV of this range, 
the DOS has a remarkably free-electron-like ^/E shape, 
without significant structure. Over the lower end of this 
region, the DOS is dominated by S 3s character, above 
which H Is and S 3p character enter and mix. Then, at -4 
eV and +5 eV two substantial and rather narrow peaks 
emerge, indicative of very strong hybridization, perhaps 
bonding and antibonding signatures. Double valleys lie 
at -2 eV and - 1-2 eV, between which a very sharp peak, 
related to two van Hove singularities (vHs) 0.25-0.30 eV 
apart in energy, juts upward. The Fermi energy Ep (set 
to zero throughout) lies very near the upper vHs. 

More details of the vHs region from WIEN2k and 
FPLO are left to Appendix A. Due to relatively small 
differences but a very sharp peak, the values of N(Ei 7 ) 
differ substantially: 0.63/eV-f.u. from Wien, 0.44/eV- 
f.u. from EPLO. Pseudopotential results will likely give 
similarly differing values. These differences arise because 
the input parameters (orbitals, sphere radii, pseudopo¬ 
tentials, cutoffs, etc.) may not be optimized for applica¬ 
tion at such reduced volumes. We demonstrate in Sec. 
IV that for physical superconducting properties, ther¬ 
mal and dynamical broadening makes details of N(E) 
fine structure relatively unimportant. This unimportance 
does not however apply for the underlying theory, where 
it has serious consequences (discussed in Sec. V) partly 
because so much is formulated and evaluated in terms of 
the specific value of N(Ei 7 ) but also because the energy 
dependence has significant impact. 

Returning to the DOS, such strong structure in N(E) 
refiects strong mixing between orbitals lying in this en¬ 
ergy range, which are the H Is and S 3p valence orbitals. 
The orbital projected DOS (PDOS) presented by Papa- 
constantopoulos et am shows that S 3s participation is 
becoming small around Ei?. Their PDOS helps to un¬ 
derstand the strong DOS structure. The peak at -4 eV 
is largely S 3p character with some H Is contribution. 
The peak at +5 eV has, surprisingly, a large contribu¬ 
tion from Bloch orbitals with d symmetry around the S 
site, with some participation of all of the orbitals besides 
S s. The peak at E^? - the important one bounded by two 
vHs - is a strong mixture of H Is with S 3p, whose cor¬ 
responding tight binding hopping parameters will have a 
correspondingly large hopping amplitude. 


In this section, a minimal tight-binding model for H^S 
will be constructed using Slater-Koster two center pa¬ 
rameters. Local basis orbitals are the atomic orbitals of 
sulfur S', Pa,, P^, Pz and the three hydrogen Is orbitals 
(s); this notation is used below. The procedure is to 
first calculate the symmetry projected Wannier functions 
(WEs), based upon the seven basis functions. The result¬ 
ing WEs at reasonably large density isocontour (shown 
in Appendix B) refiect hybridizing atomic orbitals. In 
the generation of the WEs, a set of three center hopping 
integrals is generated, with the WEs as the basis orbitals. 
Erom the WE three center integrals, simpler two center 
integrals can be obtained. Some members of this latter 
set may be over deter mined, in which case a best choice 
(average) must be made. 

The on-site energies and largest Slater-Koster param¬ 
eters are listed in Table I. Relative to Ei?=0, the on¬ 
site energies (compared to those reported by Bernstein 
et in parentheses) are: ss =-8.0 (- 8 . 6 ) eV; =-5.5 
(-5.0) eV; ep = 0.0 (-1.3) eV. The procedures used by 
Bernstein et al are not exactly the same as ours, with the 
difference indicating the level of confidence one should as¬ 
sign to these energies considering the non-uniqueness of 
tight binding representations. It is eye-catching that our 
sulfur P on-site energy is indistinguishable from E^?. 



FIG. 3: Bands from the minimal two-center tight-binding 
model described in the text (red lines), compared with the 
full band structure (black lines). The agreement is reason¬ 
able but approximate over this 30 eV region. 

The largest hopping parameters are nearest neighbor 
(n.n.) sPa (-5.4 eV) and sSa (-4.4 eV), leading to the 
possibility that H-S n.n. hopping is dominant in creat¬ 
ing the DOS structure. We have calculated the DOS for 
a single sublattice, ReOs structure H 3 S, and it is noth¬ 
ing like the full DOS. The n.n. H-H hopping ssa, at 
- 2.8 eV, couples the sublattices strongly, and S-S second 
neighbor hoppings, one lattice constant apart, are 1.3- 
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TABLE I: Slater-Koster two-center parameters (in eV).The 
subscripts indicate the index of the neighbor: l=lst neighbor, 
2=2nd neighbor, etc. For H there are two 4th neighbors one 
lattice constant apart: one through the S atom, denoted “4,” 
one through another H atom, denoted “4'.” {PP)i indicates 
|(PPcr)i + |(PP7r)i, which only occurs in this combination. 


£5 

-7.98 

{sPa)i 

-5.42 


-5.46 

{PP(t)2 

-1.83 

ep 

-0.03 

{SPa)2 

1.29 

{sSa)i 

-4.37 

{SSa)2 

0.94 

{ssa)i 

-2.80 

(sPa)2 

-0.93 

(55(7)4 

-1.14 

{PPh 

0.60 

( 55 ( 7)4 

0.55 

{SSa)i 

0.30 


1.8 eV in magnitude. Interesting are the 4th neighbor 
H-H hoppings, between atoms one lattice constant apart. 
The hopping through the S atom ( 55 ( 7)4 is -1.1 eV, while 
through an intervening H atom ( 55 ( 7)4 i^ ^i^® 

with opposite sign. The bands from this ‘minimal’ tight 
binding model are compared with the DFT bands in Fig. 

where it can be seen that it captures the general be¬ 
havior (but not the vHs) of the full band structure over 
a 30 eV range. Many more hopping parameters are nec¬ 
essary to reproduce the band structure accurately. This 
result suggests that a tight binding representation is not 
natural for H 3 S. 


C. Role of Sulfur 

What then is the role of S? We have taken another 
view of this issue, by comparing H 3 S with i.e. the S 
atom replaced by another H atom. This is simply simple 
cubic hydrogen with a lattice constant a=2.8ao, but we 
calculate it as H 3 H for comparison and uniformity. The 
DOS and band structure are presented in Appendix C. 
The DOS bears little resemblance to that of H 3 S. The 
occupied bandwidth is 15 eV, and the lower 10 eV of 
this is free electron like. A vHs is encountered at -2 eV 
followed by a remarkably linear N(E) over almost 10 eV. 
Then a second vHs signals the minimum of another free- 
electron-like high conduction band. This confirms that 
it is strong H-S sSa and sPa bonding that produces the 
strong structure in N(E) shown in Fig. 2. Thus sulfur is 
crucial in producing the DOS peak at E^? in H 3 S. We note 
that the first vHs does not appear below E^? in eitheJ^l 
fee or bee H, which are closer packed and more stable 
phases of elemental H. 

The peak in N(E) reveals that the band filling in H 3 S 
happens to be almost perfect {Ep lying very nearly at 
the sharp peak). Retaining the same band filling sug¬ 
gests substituting Se or Te for S. Flores-Livas et a/.^have 
done parallel calculations for H 3 S and H 3 Se. The H 3 Se 
frequency ujiog is 10 % higher but the calculated value of 
A is lower by 40%, with the resulting Tc being lower by 
27%. The changes of ujiog and A indicate that the product 
T] = N{Ef) < P > is lower by 20% for the Se compound. 


Here < p > is the Fermi surface average of the square 
of the electron-H ion scattering matrix element. With H 
so dominant in the EPI and H modes separated from S 
(or Se) modes, the picture is dominated by 




N{Ef) <ll> 
Mhuj% 


( 1 ) 


where the matrix element refers to scattering from the 
displaced H potential and ujh is a characteristic fre¬ 
quency from H modes. 

The other isovalent “chalcogenide” is oxygen, which is 
quite different from S chemically with H. We have found 
that the DOS of H 3 O in the H 3 S structure differs substan¬ 
tially from that of H 3 S. It may be relevant that H 2 O does 
not metalize until much higher pressures than are being 
considered here. Heil and BoerP^ have considered bond¬ 
ing, EPI, and Tc where sulfur is alloyed with other group 
VI atoms. With alloying treated in the virtual crystal 
approximation (averaging pseudopotentials), they have 
suggested that a more electronegative ion will help. This 
leaves only oxygen in that column, and they calculated 
that a strong increase in matrix elements compensates 
a considerable decrease in N(Ei?), so that A might in¬ 
crease somewhat. Ge et al. have also suggested partial 
replacemenfP^ of S, with P being the most encouraging, 
due to the increase in N{Ef); however, spectral density 
smearing will decimate this difference. 

While S changes the electronic system very substan¬ 
tially from that of H 3 H, it may not be so special. A va¬ 
riety of calculations have predicted (see Durajski et 
for references) high values of Tc (in parentheses) for H- 
rich solids: SiH 4 (H 2)2 (107 K at 250 GPa), B 2 H 6 H (147 
K at 360 GPa), Si 2 H 6 (174 K at 275 GPa), CaHe (240 
K at 150 GPa). Whether any general principles can be 
extracted from these results remains to be determined. 


D. Charge density within leV of Fermi level 

The charge density from states within 1 eV of E^? is 
shown in Fig. it is this density whose coupling to H 
vibrations gives strong coupling and the very high value 
of Tc- Results from smaller energy slices are no different, 
indicating that the states in this range have the same 
character. Two of the bands become fiat in a region 
away from symmetry lines, and the resulting two vHs 
give rise to the sharp and narrow peak. The density 
around S is strongly distorted from spherical symmetry, 
having substantial maxima in the direction of neighbor¬ 
ing H atoms. The H density is strongly elongated toward 
the two neighboring S atoms, more strongly than might 
have been guessed. These shapes reflect strong covalent 
H I5 - S 3p(j interaction, although there remains a density 
minimum in the bond center rather than a bond charge 
maximum. This character is typical of strong directional 
bonding in metallic compounds. A more nuanced indica¬ 
tion of hybridization is available from the isosurface plots 
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FIG. 4: Isosurface of the charge density obtained from states 
in the energy range Ef= 1=1 eV. Sulfur atoms (yellow spheres) 
lie at the corners and at the body center, H atom are denoted 
by small blue spheres. 


of the symmetry-projected Wannier functions presented 
in the Appendix B. 


E. van Hove singularities of H 3 S 


N(E) expanded in the energy range -1.0 eV to 0.2 
eV is plotted in Fig. 7 in Appendix B, which reveals 
the two vHs at eio and Shi from both the Wien2k and 
FPLO codes. Fitting N(E) to the following piecewise 
expressiorP^ for 3D vHs near two singularities 


N{E) 


+ Cie + di 

< a 2 e + 62 

^a3\/|e + 63I + 036 + ds 


e < Go 

^lo ^ ^ ^ ^hi 

^hi < e 


( 2 ) 


we obtain an excellent fit as can be seen in Appendix 
B, Fig. 7, demonstrating that contributions from other 
bands are smooth and slowly varying on this energy scale. 

Isosurfaces at the two vHs energies are presented in 
Fig.m The vHs points lie at either end of a line where 
low velocity regions of two sheets of Fermi surface break 
apart at the zone boundary, and then “unzip” until they 
separate into disjoint sheets. The vHs (from the FPLO 
bands) at -0.43 eV occurs at (-0.42,0.21,0)7r/a and sym¬ 
metric points. In the local principal axis coordinate sys¬ 
tem the effective masses are -O.lSrUe, 1.36me, O.ldrUe, 
giving a thermal (or DOS) mass rrith = |mim 2 m 3 |^/^ = 
0.31me. For the one at -0.11 eV, the masses are -0.83me, 
-O.lOrUe, 0.56me, and mt/i=0.42me. We return to the 
importance of vHs effective masses in Sec. VI. 



FIG. 5: Top: iso-energy contour from FPLO for the vHs at 
E=-0.43eV. Bottom: similarly for the vHs at E=-0.11 eV. The 
red circles pinpoint two of the symmetry related vHs in each 
case. The red ellipse outlines the region joined by the two 
van Hove singularities; two energy surfaces are “unzipped’ as 
the energy increases between the two vHs. Golor denotes the 
velocity, which ranges from zero (darkest tan, at the vHs) to 
2.5x10® cm/s (deep blue). 


IV. IMPACT OF ELECTRON-PHONON 
COUPLING 


Thermal broadening as described by the Fermi-Dirac 
distribution function has well understood consequences, 
and is a factor in determining (hence limiting) Tc, though 
it is usually not thermal broadening of the DOS, which is 
normally constant over an energy range of several ksTc. 
The purely thermal aspect is formalized in the gap to 
thermal broadening scale ratio 2A/7TkBTc ~ 1, and is 
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central to, but standard in, Eliashberg theory. Specifi¬ 
cally, it is thermally excited electron-hole excitations that 
overcome pairing and restores the normal state above Tc- 
However, the sharp structure in N(E) on the scale of 
relevant phonon energies requires an extension of con¬ 
ventional implementations of Eliashberg theory,^^ where 
a constant N(E) on the phonon energy scale is assumed, 
so all scattering processes can be considered as confined 
to the E ermi surface E=Ei?. This issue was confronted 
long agoj ^^ l ^^ l because of the sharp structure in N(E) in 
several of the then-high-Tc A15 structure compounds, 
viz. NbsSn, VsSi, NbsGe, with Tc^20 K, and has been 
followed up in related applications 

In H 3 S around 200 GPa, the representative frequency 
is ft 1300 K = 112 meV. We use this value below, 
based on the DET-based calculations of (harmonic) ujiog 
of 1125 K [H], 1264 K [H], 1335 K [4lT4] . and 1450 K 
[25]. The logarithmic, first, and second moments differ 
only at the 2 % level, whereas reported values differ by 
as much as 20%. Using a somewhat different value of Q 
below would not change our conclusions. 

Erom visual examination, the sheets of the constant 
energy surfaces at Ei ?±0 do not differ much from the 
Eermi surfaces at E^? (see Eig. 5), the differences occur¬ 
ring in small pockets around T (not visible) and along a 
line connecting the two vHs, irrespective of the electronic 
structure code that is used. Inter-vHs scattering could be 
interesting: although it involves a small amount of phase 
(q) space, it incorporates a disproportionate fraction of 
states with low to vanishing velocity. Possible compli¬ 
cation from inter-vHs scattering and non-adiabatic pro¬ 
cesses he beyond the scope of our discussion. 


A. Formalism 

It has been known since the work of Engelsberg 
and Schrieffei!^ and Shimojima and IchimursP^ that, 
for a characteristic phonon frequency Q in an interact¬ 
ing electron-phonon system, electron spectral density is 
spread from its noninteracting (5-function spike at Ek up 
to a few Q. The spectral density arises from the electronic 
self-energy that is treated for superconducting systems by 
Eliashberg theory. When N(E) hardly varies over a scale 
of a few Q it is rare to notice the effects of such broaden¬ 
ing except possibly in direct measurements where phonon 
sidebands may be observed in photoemission spectra. Eor 
situations as in H 3 S where Bloch state character is slowly 
varying in energy but N(E) varies rapidly, the normally 
simple electron-phonon formalism becomes challenging. 
Drozhov studied the effects of EPI in the vicinity of a 
vHs where Migdal’s theorem is violated, and found se¬ 
vere renormalizations .1^ If these are confined to a very 
small phase space, however, the effects on most proper¬ 
ties may be minor. 

Eor the case of rapidly varying N(E) but neglecting vio¬ 
lations of Migdal’s theorem, the generalization of Eliash- 
berg theory has been formulated and applied to the A15 


compoundsOne feature that is distinctive in H 3 S 
compared to most other EPI superconductors is that Tc 
is an order of magnitude higher, because the frequen¬ 
cies are comparably higher, and simple thermal broad¬ 
ening is correspondingly larger and requires attention. 
The second factor in common, and the important one, 
is that strong EPI causes an effective smearing of the 
electronic spectral density due to exchange of virtual 
phonons appearing in Migdal-Eliashberg theory - exci¬ 
tations described by the electron Green’s function are 
part electron, part phonon. This broadening is given by 
the imaginary part of the interacting electronic Green’s 
function 

= Lo-[Ek-fi{T)]-Mk{uj;T)-iTk{co-,T).{3) 

Here is the DET band energy, /i(T) is the chemical 
potential, and M and F are the real and imaginary parts 
of the phonon-induced self-energy. 

The spectral density A{uj) is the interacting analog of 
the band DOS N(E): 

= - U \ImGkiuj)\ (4) 

k k 

_ ^ _ rfc(u)) _ 

n^[LU-{Ek-fx)- Mk{uj)]^ + n(w)2 

^ j{oj-+ 

In the last expression the Brillouin sum has been con¬ 
verted into an energy integral by inserting f dES{E — 
Ek)=l, assuming that only Ek (and not wavefunction 
character, hence not M or F) depends on k near E^?, and 
ik = ^k — E M{k,^k) is the quasiparticle energy. This 
simplification is usually fine for electron-phonon coupling 
in a standard Eermi liquid, as wide-band H 3 S appears to 
be. 

There is strong rearrangement of spectral density even 
before this smearing effect of electron damping F. For 
temperatures and frequencies uj up to the order of the 
characteristic phonon energy Q or more, the behavior of 
the real part Mk is linear dMk/duo = —A/^, where Xk is the 
EPI strength at k whose average over the Fermi surface 
is A. The equation for ^k in the previous paragraph then 
gives for the quasiparticle energy 


This equation expresses the phonon-induced mass en¬ 
hancement, and (l + A/c)“^ is the quasiparticle strength, 
i.e. the fraction of the electron’s (5-function spectral den¬ 
sity at ^k and whose average in H 3 S is 1/3 (A 2 ). 
Two-thirds of the spectral weight is spread from ^k by 
up to a few times U. This is a serious redistribution of 
weight that we cannot treat in any detail without explicit 
solution for the self-energy on the real axis. 

Notwithstanding the complications, in an interacting 
system the thermal distribution function containing all 
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complexities can be handled formally to provide insight 
into this “varying N(E)” kind of system. The interacting 
thermal distribution (state occupation) function f{Ek) is 
defined as the thermal expectation of the number opera¬ 
tor rik 

^ pOO 

f{Ek) = TY,Gk{itOn)e^‘^-^ = / dujU{uj)A{Ek,co), 

-oo 

where the Matsubara sum with positive infinitesimal 77 
has been converted into an integral in the last expression 
and fo{E) is the (non-interacting) Fermi-Dirac distribu¬ 
tion. The interacting distribution function can be ex¬ 
pressed as the non-interacting one broadened b}®lr/i;(ci;) 
as N(E) is broadened in Eq. (4). 

Several thermal properties can be formulated^ in 
terms of the interacting (broadened and in principle mass 
renormalized) density of states Af{E). Returning to sin¬ 
gle particle language, the spectral density at Ef is ap¬ 
proximately 

( 6 ) 

For energies a few Q around Ef the extension Ef ^ E 
to give M{E) will be reasonable. Then, returning to the 
distribution function, the total electron number can be 
writteiP^ in two ways 

N,i = I dEf{E)N{E) = j dujfo{uj)Af{aj), (7) 

illustrating that interaction effects can be exchanged be¬ 
tween the distribution function and, in this instance, the 
interacting and non-interacting density of states. Thus 
in a region around Ef the spectral density is the band 
density of states broadened by a Lorentzian of halfwidth 
T. 


B. Thermal and phonon smearing in H 3 S 

We now estimate the impact of this spectral density 
smearing for H 3 S using the Wien 2 k result for N(E). The 
mass renormalization effects (from the real part of the 
self-energy) Migdal theory are outlined in Appendix D 
but will be disregarded here, leaving our estimate as an 
underestimate of the effect of smearing. Investigation 
of the (Migdal) self-energy equations (ie. in the nor¬ 
mal state) gives the quasiparticle inverse decay r ate via 
phonons over most of the relevant energy range,^ ^^ * ^^ * 
for an Einstein model, as 

r PS TT\fi[nB{ei) + -] (8) 

For H 3 S the characteristic frequency Cl is in the ballpark 
of Cl ^ 1300 K (see above). Since we are interested only 
in relative values of quantities affected by smearing, we 
will not distinguish ujiog from UJ 2 ^ etc. 


With the choice of /i*=0.15, a value of A=2.17 is nec¬ 
essary to give the observed Tc=200K, which also is in 
the range that has been quoted as resulting from DFT 
calculations. At 200 K, the Bose-Einstein thermal distri¬ 
bution ub gives a negligible fraction of phonons excited: 
nB{Cl) - 10-^. Thus 

r = = 5 X 10^ X = 0.38 eV, (9) 

the halfwidth is proportional to the product XCl. This 
smearing in Eq. (6) arises from the virtual excitation 
of phonons that provides the coupling, even at low Tc 
where phonons are not excited. We note that zero-point 
vibrations do not scatter electrons; sharp Fermi surfaces 
survive strong electron-phonon coupling, with de Haas - 
van Alphen oscillations remaining visible, and the resis¬ 
tivity p{T) ^ 0 at T ^ 0. 



Energy (eV> 

FIG. 6 : H 3 S density of states without broadening (black line, 
Wien 2 k DOS), with thermal broadening at 200 K (green line, 
which is hard to distinguish from the unbroadened one), and 
the virtual phonon broadened effective DOS N(E) using a 
Lorentzian halfwidth of 0.38 eV (red line). Note the large 
drop in the value at the Fermi energy (dashed vertical line). 

Figure]^ shows three curves: the calculated (Wien 2 k) 
static lattice N(E) from Fig. 2, the thermally broad¬ 
ened version at 200 K, and the virtual phonon broadened 
DOS with halfwidth F=0.38 eV from above. The thermal 
broadening function is the derivative —dfo{E — E')/dE 
with half width of about Trk^T = 625 K = 55 meV. This 
thermal broadening at 200 K is minor on the scale of 
interest, reducing the effective value of N(Ei 7 ) slightly. 
The phonon broadening however is severe, with the peak 
value of 0.70/eV-f.u. for N(E) dropping by 37%. The un¬ 
broadened value N{Ef) of 0.64 states/eV-f.u. is lowered 
by nearly 1/3, to 0.45 states/eV-f.u. The effect of the 
shift of the chemical potential is secondary when broad¬ 
ening is so large. 
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C. Implications for the theory of H 3 S 

The El iashbe rg equations including the energy 
dependenc d^^ l ^^ l of N(E) indicate that it is this 1/3 re¬ 
duced value of N(Ei 7 ) that should be used with the 
standard implementation to get a good estimate of A, 
the superconducting gap, and Tc- The impact of ex¬ 
periment/theory agreement for H 3 S is substantial and 
negative: the spectral function a^F(co’), proportional to 
N{Ef)^ is reduced by 1/3 by EPL The naive value of 
A ~ 2.17 becomes, after reducing by 1/3, A=1.45. The 
phonon frequency moments, which involve a^F/A, seem 
at this stage to remain unchanged; we return to this point 
below. 

We have evaluated the magnitude of this phonon 
broadening correction on Tc using the Alien-Dynes 
equation,!^ taking the representative values for H 3 S of 
the phonon moments to be (1 = 1300 K and yU*=0.15. Eor 
A=2.17, Tc = 200 K; for the 1/3 reduced value A=1.45, 
Tc=130 K. The agreement between theory and experi¬ 
ment is strongly degraded. It is worthy of note that H 3 S 
happens to be in a nearly linear regime of Tc(A), where a 
reduction of A by 1/3 results in a decrease of Tc by nearly 
1/3. We note that the corresponding strong coupling fac¬ 
tor fi in the Alien-Dynes parametrization is 1.13 (1.07); 
fi is the crucial improvement of the Alien-Dynes equa¬ 
tion over the McMillan equation and seems sometimes 
for hydrogen sulfides to have been neglected. 

What this comparison implies is that theory- 
experiment agreement for Tc is not as good as has 
seemed, since taking phonon smearing into account, the¬ 
ory would only be predicting of the order of 2/3 of the 
“constant A^(Fi 7 )” value. It should be noted that this 
1/3 reduction factor depends on the accurate calculation 
of N{E)^ for which we have used the Wien2k spectrum. 
With the EPLO result for the DOS, N{Ef) is lower and 
thus the effect of smearing will be smaller. 

A few papers have reported calculations of the band 
structure and of N(E), and some studies have reported A, 
but little attention has been given to the value of N(Ei?), 
which is sensitive to method and computational proce¬ 
dures. At 200 GPa (a=5.6 a.u.) Papaconstantopou- 
los et al. quote 0.51/eV-f.u.; from the figure of Duan 
et al. we estimate the small value of 0 . 2 /eV-f.u. (cer¬ 
tainly their presented peak is much weaker): Bianconi 
and Jarlborg report 0.50/eV-f.u. (their table numbers 
must be per cubic cell). Our Wien 2 k and EPLO values 
are 0.64 (0.42)/eV-f.u. respectively, indicating that even 
all-electron full potential methods can differ. 

The point is that in conventional Eliashberg theory - 
constant N(E) on the phonon scale - A is proportional 
to N(Ei 7 ), and the values that have been used are sensi¬ 
tive to methods and cutoffs (depending on method, see 
above), but more seriously they are obtained from un¬ 
broadened N(E). Because of this, the reported values of 
A and hence Tc are quantitatively uncertain, assuming 
they are converged BZ integrals. And on this point, 
McMahon and Ceperle^ttSl and Akashi et have dis¬ 


cussed the various challenges in reaching convergence, be¬ 
fore even confronting the energy variation question. The 
(unsmeared) prediction of Tc~ 200 K indicates that im¬ 
proved theory, by taking into account phonon broaden¬ 
ing, would give a substantially reduced critical tempera¬ 
ture. 

Elores-Livas et al have recognized the issue of the vari¬ 
ation of N(E) on the scale of the phonon frequencies. In 
their implementation of density functional theory for su¬ 
perconductors (DETSC), this variation is accounted for. 
They reported, for their calculation of N(E) (details were 
not reported), taking into account the variation resulted 
in a 16% decrease in Tc, from 338 K to 284 K in their 
calculation. Their methods also involve calculation of /i* 
that is not treated in the Eliashberg form as well as other 
methodological differences, and these differences make di¬ 
rect comparison with other reports difficult. Still, the 
relative effects of energy variation of N(E) are clear. 

This DOS variation issue extends to the calculation 
of phonon frequencies. The phonon self-energy involves 
electron-phonon scattering in which a phonon is ab¬ 
sorbed, scattering an electron from Ef^ < E^? to E}^-\-ujq > 
Ef- Most methods of calculating phonon frequencies do 
not include effects of the density of available initial or fi¬ 
nal states in this energy region being variable. Thus cal¬ 
culation of phonon spectra will need to be re-evaluated 
for situations such as that imposed by H 3 S, and the as¬ 
sociated non-adiabatic corrections considered. 


V. SCENARIOS FOR ROOM TEMPERATURE 


The foregoing section indicates that 200 K supercon¬ 
ductivity has been achieved with an effective DOS around 
M{Ef) ~ 0.44/eV-f.u. compared to a peak value be¬ 
tween the two vHs around 0.7/eV-f.u. This means as 
mentioned above that the theory needs refining, as al¬ 
ready noted by Elores-Livas et a/., to determine just how 
much is understood quantitatively and what features may 
require more attention. 

This difference in N(E) versus M(E) has larger and 
much more positive implications. The fact that a re¬ 
duced effective value of M{Ef) should be used for H 3 S 
also provides the glass-half-full viewpoint: a much larger 
value of M{Ef) and therefore Tc may be achievable in 
this or similar systems. Suppose that the two vHs can be 
moved apart, each by (say) 0.5 eV, leaving a value N(Ei 7 ) 
~0.7/eV-f.u. between. Then the Lorentzian smearing 
will have much less effect. There is also the question 
of increasing the magnitude of N(E) at the peak, i.e. 
N{Eyhs)- This value depends on the effective masses at 
the vHs, but also on the volume of the region in which 
the quadratic effective mass representation holds. If it 
holds in an ellipsoidal region defined by 
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the DOS from this region is 

+ ...], ( 11 ) 

where a and P are numbers of order unity. The value 
N(E^/is) at the vHs is proportional to the thermal mass 
and to the radius Gc of the region of quadratic disper¬ 
sion, and the decrease away from the vHs (the second 
term) is inversely proportional to Gc- There will be ad¬ 
ditional smooth contributions from outside this region, 
of course. However, increasing rrith and the region of 
quadratic dispersion is favorable for increasing Af{E) in 
the vHs region, and hence increasing Tc. These ob¬ 
servations seems to implicate the topology of the Fermi 
surfaces, rather than more conventional electronic struc¬ 
ture characteristics such as relative site energies and hy¬ 
bridization strengths. 

Numerical examples are illuminating. Suppose that 
the DOS peak can be widened so that M{Ef) N(Ei?) 
= 0.70/ev-f.u. as outlined just above, rather than the 
reduced effective value of 0.45eV/f.u. that gives, exper¬ 
imentally, Tc=200K. With ft = 1300 K and /i*=0.15 as 
above, Aea:p=2.17 is required to account for Tc = 200K. 
For a 0.70/0.44 = 1.55 larger value of N(Ei 7 ), A=3.38 and 
we find Tc=277K - room temperature in a cool room. 
The increase in effective M{Ef) we have assumed is am¬ 
bitious but not outlandish, given the calculated spectrum 
of ImSm H 3 S. It is clearly worthwhile to explore other 
H-rich compounds for higher critical temperatures. 

Of course, the increase in JV{Ef) will give additional 
renormalization (softening) of the phonons. However, the 
modes are very stiff even with A= 2 . 1 , so this may not be 
a major effect. Note that decreasing ft increases A but 
decreases the energy scale prefactor in Tc, one reason 
why increasing A by decreasing frequencies is rarely a 
profitable means of increasing Tc. If fl, the prefactor in 
the Alien-Dynes Tc equation, is softened by 10% with¬ 
out change in matrix elements, A increases by 20% while 
Tc increases by only 4%. Evidently softening of hard 
phonons is a minor issue when looking for higher Tc in 
this range of A. This behavior was formalized by Allen 
and Dynes ,1^ who obtained the rigorous strong coupling 
limit 

Tc O.ISVA < > = 0.l8^y N{Ef) <P> /M, (12) 

where uj is expressed in keIvins. The last expression is 
strictly true only for an elemental superconductor, think 
of the electron-ion matrix element / and mass M as those 
of H for H 3 S. 


VI. SUMMARY 

In this paper we first addressed the electronic structure 
and especially the delicate van Hove singularity induced 
spectrum, bonding characteristics and effect of S, and 
the charge density of states near the Fermi level, more 


directly than has been done before. The occurrence of 
two closely spaced van Hove singularities is definitely a 
central issue for the properties of H 3 S. We list some of 
the main points. 

• At the most basic level, why is H 3 S superconduct¬ 
ing at 200 K? It is because both A is large but, more 
importantly, the characteristic phonon frequency ft is 
very high. This reminds one of the Alien-Dynes limit 
for strong coupling, 

Tc 0.18\/A1F ^ 0.18y/N{EF) <P> /M. (13) 

Though not yet in this limit, this provides the right pic¬ 
ture - one can check that keeping all fixed except for ft 
and then varying it, the change in Tc is minor because the 
change in prefactor Tc (x ft is compensated by A (x ft~^. 

• Sulfur 3p states hybridizing with hydrogen Is is cru¬ 
cial in producing the strong large scale structure in N(E) 
within ±5 eV of the Fermi level, and in leaving Ep at the 
top of a peculiarly sharp peak between two vHs. The van 
Hove points on the constant energy surfaces that define 
the peak in N(E) were identified, finding they lie on op¬ 
posite ends of a line of Fermi surface “ripping apart” with 
energy varying between the two van Hove singularities. 
This region of very low velocity electrons affects a sig¬ 
nificant fraction of the zone. It is unclear how replacing 
S with other elements will affect the electronic structure 
near E^?, but small changes may have large effects. Ge 
et al. have noted that alloying 7-10% of P with S moves 
Ep to the peak in N(E), within the virtual crystal ap¬ 
proximation which does not account for alloy disorder 
broadening. In any case, strong coupling smearing as 
discussed here will nullify this apparent gain. 

• The fine structure and energy variation of N{E) near 

the Fermi level must be taken into account to obtain 
quantitative results for A, and Tc. The energy de¬ 

pendence of N(E) may even affect calculation of phonon 
frequencies, though this is untested so far. 

• The closely spaced van Hove singularities very near 
the Fermi level have been shown to have significance, 
both on the detailed theory of H 3 S but, as importantly, 
on the question of whether Tc can be increased in related 
materials. Sulfur and the specific ImSm structure are 
important for high Tc though other elements will need to 
be studied to learn more about precisely why. 

• The prospect for increased Tc is affirmative - it 
will require only evolutionary changes of the electronic 
structure to achieve room temperature superconductiv¬ 
ity, though the road to this goal requires study, and ad¬ 
ditional insight into the origins of van Hove singularities 
may be important. Increasing the vHs effective masses, 
or increasing the volume within which quadratic disper¬ 
sion holds, will increase N(E) at the vHs energy. Struc¬ 
tural or chemical changes that affect the electronic struc¬ 
ture rather modestly may lead to significant increase in 
the effective (broadened) density of states at Ei?. Other 
studies have suggested that substitution of some sulfur 
with chemically related elements may increase Tc. Alto- 
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gether, the prospects of achieving increased critical tem¬ 
peratures are encouraging. 

• An issue that is almost untouched is a deeper under¬ 
standing, or rather an understanding at all, of electron- 
ion matrix elements < Ijj > - what contributes to strong 
electron-H atom scattering, and what degrades this scat¬ 
tering. These matrix elements are the same that de¬ 
termine resistivity in the normal state; notably most of 
the best superconductors have high resistivities. Further 
study should address the EPI matrix elements. 
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Appendix A: Details of van Hove singularities 


In FPLO, which we used initially to identify the effec¬ 
tive masses at each vHs, the vHs energies he at -0.43 eV 
and -0.11 eV, i.e. the Fermi level lies 110 meV above 
the upper vHs. In Wien2k, the vHs lie at -0.20 eV and 
+0.05 eV; Ei? lies just below the upper vHs, thus within 
the peak. The values of N{Eyhs) ^.t the upper vHs are 
0.555/eV-f.u. (FPLO) and 0.630/eV-f.u. (Wien2k), re¬ 
flecting the difference in thermal masses. Pseudopoten¬ 
tial results may differ by somewhat more than do these 
two methods, which are usually in excellent agreement. 


Appendix B: Minimal tight-binding model 


A minimal tight-binding model for H^S has been con¬ 
structed using Slater-Koster two-center hopping parame¬ 
ters. To simplify notation, we denote the sulfur 3s orbital 
as S', 3p orbitals as P, and hydrogen Is orbital as s. The 
subscripts indicate the neighbor of the second site rela¬ 
tive to the first. The Slater-Koster parameters, provided 
in Table I, are discussed in the text. 



FIG. 7: Expanded view of the van Hove singularity region 
for WIEN2k (top panel) and FPLO (bottom panel) calcula¬ 
tions. In this energy range, there are two prominent van- 
Hove-singularities within 0.25 to 0.30 eV of one another, with 
the top one lying very near Ef- The red dashed lines mark 
the Fermi level. Red, blue and green dashed lines are the fit¬ 
ted DOS using Eq. N{Ef) from FPLO is 0.44 states/eV, 
while from WIEN2k it is 0.64 states/eV. 


indices, 1 is S S, 2-4 are S P, and 5-7 are the H s orbitals. 


Hu = {SSa)o E‘2{SSa)2{coskxcosky + cosk(^T^l) 

Hi2 = 2i{SPa)2smkx 


2{Ssa)i cos ^ 


Pl5 

H 22 = 2(PP)i[cos 


kx + ky + k^ 


kx + ky 


ky -\- kz 


+- —] + 2 {PP(j )2 cos kx 


/ -I—« \ / / . kx kz kx I kz \ 

H 26 = -{sPa) 2 V 2 t{sm - - -h sin- - -) 

^55 = (sscr)o + 2(sSCr)4 COS/Ca; + 

2(sscr)4 COS ky + 2{ssa4)' cos kz 

P56 = 2 (sSCr)iCOSy 


Sulfur Px to Py hopping is small to the nearest neigh- 


Selected S-K matrix elements are provided that, by 
permutation, will allow construction of the tight binding 
model in the form we have used. For the row and column 
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FIG. 8 : Left: Sulfur 3pz Wannier function, showing strong 
mixing with H Is states above and below. Center: Sulfur 3s 
Wannier function, with minor hybridization to neighboring Is 
orbitals. Right: Hydrogen centered Wannier function, reveal¬ 
ing strong hybridization of H Is with S 3p above and below,, 
as well as mixing with neighboring H Is orbitals along the 
other two axes. 


bor and vanishes by symmetry for the second neighbor, 
so H23=0 at this level of the model. The dispersion ob¬ 
tained from the model Hamiltonian gives a reasonable, 
but not quantitatively accurate, representation of the 
DFT bands. 


Appendix C: H 3 H electronic structure 

The density of states and band structure of the model 
compound H 3 H discussed in the text are presented here. 
This “compound” is actually simple cubic hydrogen, but 
presented in the H 3 S type cell to facilitate compari¬ 
son. The occupied bandwidth is 15 eV. In bee and fee 
structures,^^ there is no vHs below Ep. The conclusion 
is that sulfur has a momentous impact on the electronic 
structure within 5 eV of the Fermi level. 



FIG. 9: H 3 H density of states using WIEN2k with Rh = 
1.39ao. RnKmax is set to 6 to be consistent with the H 3 S 
calculations. The k-mesh was 100x100x100 with 22750 k- 
points in the irreducible wedge, and the tetrahedron method 
was used in generating the density of states. 


[[3U W[E'N2kl’iandsLmtlurc 



FIG. 10: H 3 H band structure from WIEN2k, using the cell of 
H 3 S and calculational parameters as in the caption to Eig. 


Then returning to Eq. (4) before neglecting M/^, 


Appendix D: Spectral redistribution including Mk{oj) 


Here the lower frequency region of the smearing of 
the spectral density is outlined. We consider the regime 
where dMkiuo)/duj = — A/^, which holds up to ~ though 
perhaps not as far as F = ( 7 r/ 2 )All, before the renor¬ 
malization (mass renormalization) begins to “burn off.” 


= Er— 


rfc(w)/7 


[w - {Ek - jj) - Mfe(w)]2 + Tkito) 






rfc(w)/7 


[{l + Xk)co-{Ek-fiW+Tk{uj)^ 




TfcMMl + A) 


(w- 6)]2 + rfeH2/(i + A) 


(1 + A)-3/2| cimo 


f/TT 


{to - C )2 + r 2 
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The identity 1 = J — ^k)d^ has been inserted, = 
(^/c —/i)/(l + A/c), and in the last line /^-dependence other 
that through has been averaged. 

In this expression the width T = T/ y/l + A has been 
decreased by mass renormalization by a factor of 1.8 for 
A=2-2.2, reducing smearing. However, the overall mag¬ 
nitude has been decreased by (1 + A)“^/^ ^ 1/6, a very 
large expulsion of spectral weight from low energy by 


the mass enhancement. This expression provides addi¬ 
tional insight into the substantial redistribution of low 
energy spectral weight beyond the Lorentzian broaden¬ 
ing (in the text) that takes over for uj > Q. Calculations 
of the self-energy for sulfur hydride^^ show that M{uj) 
differs only moderately from its low frequency form for 
uj up to roughly the 2 VL range. 
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